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ABSTRACT 

We have determined the theoretical predictions for the cross-sections of squark and gluino 
production at pp and pp colliders (Tevatron and LHC) in next-to-leading order of su- 
persymmetric QCD. By reducing the dependence on the renormalization/factorization 
scale considerably, the theoretically predicted values for the cross-sections are much more 
stable if these higher-order corrections are implemented. If squarks and gluinos are dis- 
covered, this improved stability translates into a reduced error on the masses, as ex- 
tracted experimentally from the size of the production cross-sections. The cross-sections 
increase significantly if the next-to-leading order corrections are included at a renormal- 
ization/factorization scale near the average mass of the produced massive particles. This 
rise results in improved lower bounds on squark and gluino masses. By contrast, the 
shape of the transverse-momentum and rapidity distributions remains nearly unchanged 
when the next-to- leading order corrections are included. 
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1 Introduction 



The supersymmetric extension of the Standard Model Q is a well-motivated step. In 
super symmetric theories the hierarchy problem of the Higgs sector can be solved @. 
Even in the context of very high energy scales, as required by grand unification, it is 
possible to retain fundamental scalar Higgs particles with low masses. This is a conse- 
quence of pairing bosons with fermions in supersymmetric multiplets, which removes the 
quadratic divergences due to these high scales from the quantum fluctuations. Moreover, 
the electroweak Higgs mechanism can be generated radiatively ||. For a top mass in 
the experimentally observed range, the theory can evolve from a symmetric phase at the 
grand-unification scale to a phase of broken electroweak symmetry at low energies, while 
leaving the electromagnetic and colour gauge symmetry unbroken. Strong supporting 
evidence for supersymmetry is provided by the successful theoretical prediction of the 
electroweak mixing angle @], sin 2 9 W = 0.2334 ±0.0035, based on the particle spectrum of 
the Minimal Supersymmetric extension of the Standard Model (MSSM). This prediction 
is matched quite well by the measured value ||, sin 2 6>° xp = 0.2317 ± 0.0004. Supersym- 
metric extensions offer solutions for many other problems that cannot be solved within 
the Standard Model [see e.g. Ref. |J for a comprehensive discussion]. 

In the MSSM [0] quarks and leptons are paired with squarks and sleptons, gauge and 
Higgs particles with gauginos and higgsinos. Supersymmetric QCD (SUSY-QCD) is based 
on the coloured particles of this spectrum: quarks and spin-0 squarks (q = (Jl, qn), gluons 
and spin-1/2 gluinos (g). The magnitude of the SUSY-QCD interactions is set by the 
gauge and Yukawa couplings g s and g s = g s , respectively; the two couplings are required 
by supersymmetry to be equal. If supersymmetry were an exact symmetry, squarks and 
quarks would have equal masses and gluinos would be massless. However, supersymmetry 
is a broken symmetry, and the masses of the supersymmetric partners must exceed the 
masses of the Standard-Model particles considerably. Even though the mechanism for 
the breaking of supersymmetry is not identified yet at the fundamental level, requiring 
that no quadratic divergences be reintroduced into the theory by breaking the symmetry, 
provides a powerful guiding principle |§. In this approach, supersymmetry is broken 
within SUSY-QCD by introducing heavy masses for squarks and gluinos, lifting the mass 
degeneracy with quarks and gluons. In order not to ruin the solution of the hierarchy 
problem, the masses of these particles should not exceed limits of 0(1 TeV). 

In the present analysis we will assume that the scalar partners q = {qLiQR) of the 
five light quark flavours are mass degenerate. We defer the discussion of final-state stop 
particles, with potentially significant L R mixing due to the large top-Higgs Yukawa 
coupling, to a subsequent paper |§. Stop masses in loop diagrams can be identified with 
the other squark masses since their impact is small. The masses of the five light quarks are 
neglected, while the top-quark mass is set to m t = 175 GeV. [The set of Feynman rules 
in SUSY-QCD, including an extensive discussion of the Majorana gluinos, is summarized 
in Appendix 0; more details are given in Ref. [ 10 1 .] 



The search for supersymmetric particles ranks among the most important experimen- 
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tal endeavours of high-energy physics. The coloured particles, squarks and gluinos, can 
be searched for most efficiently at hadron colliders. As R parity is conserved in the QCD 
sector of N = 1 supersymmetric theories, these particles are always produced in pairs. 
Squarks and gluinos decay primarily into cascades of jets plus the lightest supersymmetric 
particle (LSP), which escapes undetected. This particle is in general identified with the 
lightest neutralino Xv Missing momentum is therefore one of the classical characteris- 
tics of SUSY events. Pairs of gluinos can decay into like-sign dileptons plus two LSPs, 
providing an almost background-free signal. 

Squarks and gluinos can at present be searched for at the Fermilab Tevatron, a pp 
collider with a centre-of-mass energy of 1.8 TeV, which will be upgraded to 2 TeV in 
the near future. Lower bounds on the squark and gluino masses have been set by both 
Tevatron experiments, CDF and DO [12| . At the 95% confidence level, the lower 
bound for the gluino mass is 173 GeV, independent of the value of the squark mass. If 
squarks and gluinos have the same mass, the lower limit is given as 225 GeV. Within 
the theoretical set-up of the experimental analysis, no limit for squark masses can be 
derived if the gluino mass exceeds 550 GeV. [The lower limits on squark masses obtained 
at LEP are independent of this assumption.] This set of experimental bounds on squark 
and gluino production has been obtained in the framework of supergravity models, in 
which gluinos cannot be much heavier than squarks. If the supergravity relations are 
relaxed to the supersymmetric GUT relations between the gaugino masses, the excluded 
range in the gluino/squark-mass plane can be extended slightly. For gluino masses above 
550 GeV, no bound on the squark masses can be obtained any more, since the LSP 
Xi becomes so heavy that the missing transverse momentum is insufficient to generate 
an observable signal. [Very light gluinos, which may have escaped detection at hadron 
colliders, are improbable in the light of the observed topologies of hadronic Z decays; see 
e.g. Ref. fH|]. In the near future the search for squarks and gluinos can be extended at 



the upgraded Tevatron to masses between 300 GeV and 400 GeV. At the pp collider LHC 
the mass range up to 1.5-2 TeV can be sweeped, which covers the canonical range of the 
coloured supersymmetric particles, yet may not exhaust the parameter space entirely. 

The cross-sections for the production of squarks and gluinos in hadron collisions were 
calculated at the Born level already quite some time ago ||14|1 . Only recently have the 
predictions been improved by next-to-leading order (NLO) SUSY-QCD corrections for 
squark-antisquark [15| and gluino-pair production Jl6| in pp collisions^. In the present 



paper we give a systematic analysis of the next-to-leading order SUSY-QCD corrections 
of all possible supersymmetric pair channels, 

pp/pp — > qq, qq, qg, gg + X (q^t) (1) 

in the proton-antiproton collisions at the Tevatron and the proton-proton collisions at 
the LHC. 



2 The gluonic radiative corrections to gluino-pair production in gluon fusion are closely related to the 
gluonic corrections for heavy-quark production [fj"7| , just requiring the appropriate change of the Casimir 
invariants. 
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Several arguments demand the NLO SUSY-QCD analysis in order to obtain adequate 
theoretical predictions for the cross-sections: 

(i) The lowest-order cross-sections depend strongly on the a priori unknown renormal- 
ization/factorization scale. As a result, the theoretical predictions are uncertain within 
factors of 2. By implementing the NLO corrections, this scale dependence is expected to 
be reduced significantly. 

(ii) Drawing from the experience with similar hadronic processes [e.g. hadroproduction 
of top quarks] , the NLO corrections are expected to be positive and large, thus enhancing 
the production cross-sections and raising the presently available (conservative) bounds on 
squark and gluino masses. 

(Hi) When squarks and gluinos will be discovered, the comparison of the measured total 
cross-sections with the theoretical predictions can be used to determine the masses of 
the particles. Due to the two escaping LSPs that are produced in the final-state decay 
cascades, the masses of squarks and gluinos cannot be determined by reconstructing 
the original squark and gluino states. [Transverse-momentum spectra can eventually be 
exploited to evaluate squark and gluino masses from the final-state distributions JIB] ]. 

Because of the large number of mechanisms, the calculation of the NLO corrections is 
tedious but straightforward. The only non-straightforward component of the theoretical 
set-up is one element of the renormalization program. To make maximal use of the 
common infrastructure developed earlier for top-quark production ]17j, we have chosen 



the MS renormalization scheme. However, this scheme leads in n ^ 4 dimensions to a 
mismatch between the 2 fermionic gluino degrees of freedom and the (n — 2) transverse 
gluon degrees of freedom. As a consequence, the gauge coupling g s (MS) and the Yukawa 
coupling g s (MS) of the MS scheme differ in higher orders by a finite amount, even in 
exact supersymmetric theories. The problem, however, can be solved by introducing a 
proper counter term that restores the supersymmetry also in higher orders [|T9|j . 



The paper is organized as follows. In the next section we recapitulate the lowest-order 
cross-sections for the partonic subprocesses of squark and gluino production for the sake 
of completeness and to define the notation. In Section 3 we carry out the calculation of 
the virtual SUSY-QCD corrections, followed by real-gluon radiation and the discussion of 
final states including light quarks. In Section 4 we present the overall corrections at the 
parton level and at the hadronic level for the total pp and pp cross-sections. Furthermore, 
the transverse-momentum and rapidity distributions for semi-inclusive squark/gluino final 
states will be discussed briefly. We conclude the paper with an assessment of the results. 
Some useful technical details are collected in several Appendices. 
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2 Squark and Gluino Production in Leading Order 



qq production: 


qi+qj - 


-> qk+qi 




9+9 - 


-» ?< +qi 


qq production: 


Qi+Qj ~ 


— > cji +qj and c.c. 


gg production: 


qi+qi - 


^9+9 




9+9 - 


^9 +9 


qg production: 


qi+g - 


— > q~i +g and c.c. 



To set the stage for the subsequent discussion of higher-order effects, it is useful to reca- 
pitulate the lowest-order processes of squark and gluino production in quark and gluon 
collisions Hl4| . Moreover, the main features of the production mechanisms will be briefly 
described. 

The hadroproduction of squarks and gluinos in leading order (LO) of the perturbative 
expansion proceeds through the following partonic reactions: 

(2) 
(3) 
(4) 
(5) 
(6) 
(7) 

The momenta of the two partons in the initial states are denoted by k\ and those of 
the particles in the final states by pi and p 2 . 

The chiralities of the squarks q = (q~L, Qr) are not noted explicitly. The indices i—l 
indicate the flavours of the quarks and squarks. Also charge-conjugated processes (c.c.) 
are possible, related to the reactions (4) and (7); for the sake of simplicity they are 
not given explicitly. They will be taken into account properly when the hadronic cross- 
sections are calculated. The Feynman diagrams corresponding to these partonic reactions 
are displayed in Fig. |l| The production of squark-antisquark final states requires quark- 
antiquark (a) or gluon-gluon (b) initial states. Squark pairs can only be produced from 
quark-pair (c) initial states. Gluino pairs are produced from quark-antiquark (d) and 
gluon-gluon (e) initial states. The squark-gluino final states can only be produced in 
quark-gluon (f) collisions. 

We use the Feynman gauge for the internal gluon propagators. For the external gluons 
only the transverse polarization states are generated. In the axial gauge the polarization 
sum for the external gluons is given by 



j, [rii hi) [rii kiY 



Here n« ^ hi is an arbitrary vector. This polarization sum obeys the transversality 
relations 

L. pt*v _ p^U- —n- P^—P^n- — D 
"'t/i 1 i 1 % lb i[L 1 i 1 % 'His u - 
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Figure 1: Feynman diagrams for the production of squarks and gluinos in lowest order. 
The diagrams without and with crossed final-state lines [e.g. in (b)] represent t- and u- 
channel diagrams, respectively. The diagrams in (c) and the last diagram in (d) are a 
result of the Majorana nature of gluinos. Note that some of the above diagrams contribute 
only for specific flavours and chiralities of the squarks. 
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Combining these transversality relations with the Slavnov-Taylor identities for on-shell 
external particles, e.g. 



k?M, 



MP 



Mf ost oc k jp 



(9) 



for two external gluons [labelled i and j] , allows us to perform ghost subtraction. Because 
of the transversality relations all terms proportional to fcf and kj can be removed from 
the LO matrix elements, resulting in the nullification of the right-hand side of Eq. (|9]) 
[ghost subtraction]. As a consequence, the polarization sum can effectively be replaced 

by PI 



fJIV 



The rif = 5 light quark flavours and the gluons are treated as massless particles. Since 
we have excluded the top-squarks from the final state, all squark-flavour and chirality 
states are considered to be mass degenerate with mass rriq. The gluino mass is denoted 
by rrig. The set of kinematical invariants used for the description of the reactions (2)— (7) 
is given by 



s = (ki + k 2 ) 2 
t = (k 2 -p 2 ) 2 
u = (hi -p 2 ) 2 



(10) 



ti = (k 2 - p 2 ) 2 - m\ 
ui = {fa -p 2 ) 2 - mj 



{k 2 - p 2 f - m 
{fa -p 2 ) 2 - m 



The Mandelstam invariants are related by s+t+u = p\+p\. All in- and outgoing particles 
are assumed to be on their respective mass shells, i.e. k 2 = 0, p 2 = m~ for squarks, and 
Pi = m 1 f° r ghiinos. 

Applying the Feynman rules given in Appendix we obtain the following squared 
lowest-order matrix elements l-M' 8 ] 2 in n = 4 — 2e dimensions^: 



Y,\M B \ 2 {m^m) 



8n f g 4 NC F 



t\Ul 



-^- + Ag 4 s NC F 



(m~ — m~)s 



8g 2 g 2 c F - 



m 2 ~s 



si,, 



+ (1 - S i: 



Ag 4 s NC F 



3 

hui - (mj - mps 
t 2 



*£\M B \ 2 {gg^fq) 



Co 1 



■Ok 




sm 2 - 



Y,\M B \ 2 { qi q 3 ^qq) 



2g 4 s NC F {t lUl 



m 2 gs) 



1 1 

To o 



t 2 



II; 



+ I/A mi* [ AYV (^ + ^]-2C F ' 



+ (1 - 6. 



Ag 4 s NC F 



t 2 u 2 
hut - {m 2 , - m])s 



tgUg 



3 The generalization to n ^ 4 dimensions is anticipated here in view of the renormalization and mass 
factorization that will have to be performed in higher orders later on. 
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"£\M B \ 2 (qq^gg) = lg A a C 



^2\M B \ 2 (gg^gg) 

J2\M B \ 2 (gg^q~g) 



2m 2 ~s + t 2 g + u 2 g 



+ Ag 2 s g 2 s Co 



m 2 g s + t 2 g | m 2 g s + u 2 g | J tg | u g 

St\ SU\ \t\ U\ 



+ 2gt 



Co 



it 



+ 



+ C K 2 



mis 



t\ U\ ti 



SgtNColl-^ 



s 2 9 fn-s 

(l-e) 2 -2(l-e)+4 ; s 



tg Ug 



tgUg 



0,|J-2^ 



<- 1 + E >7 



2(m| - m|) t g 



1 + 



m 2 g s 

tgUg 



The QCD gauge coupling (ggg) is denoted by g s and the Yukawa coupling (qqg) by g s . 
These couplings are identical. For squarks all chiralities and non-stop flavours are summed 
oveif]. As mentioned before, the charge-conjugate final states that are not denoted ex- 
plicitly in the above listing, will be included in the hadronic cross-sections. The SU(3) 
colour factors are given by iV = 3, C = N(N 2 - 1) = 24, C K = (N 2 - 1)/N = 8/3, and 
C F = (N 2 - 1)/(2N) = 4/3. 

After performing the n-dimensional phase-space integration and taking into account 
colour and spin averaging, we find for the lowest-order double-differential distributions: 



S dtdu ' n r(i -, 



(t p 2 ) p\s 2 2 2 

®([t ~ P2\[ U ~ P2\ -P2S 



jJ?S 



x 0(a - 4m 2 ) 5(s + t + u - p\ - p 2 2 ) ^ \M B \ 2 . (11) 

Here m denotes the average mass of the produced particles, i.e. m = (\fp{ + \/p§)/2. 
The averaging of the initial-state colours and spins is incorporated in the factor Kif 

K = K - = K = - K ' 



4iV 2 ' 99 4(l-e) 2 (iV 2 -l) 2 ' q9 4(1 - e)iV(iV 2 - 1) ' 

The gluons have (n — 2) degrees of freedom in n dimensions. The scale parameter /x 
accounts for the correct dimension of the coupling in n dimensions. The term S £ = 
(47r) _2+e follows from the angular integrations. 



4 The first term of YJ \M. B \ 2 {q%qj — + qq) corresponds to the production of squarks with different 
chiralities, whereas the second term corresponds to equal chiralities. When calculating cross-sections, 
the second term will be weighted by a statistical factor of 1/2, since the squarks in the final state are 
identical. 
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The subsequent integration over the remaining invariants^] yields the total lowest-order 



partonic cross-sections |L4 

v B (<nqj -> qq) : 



° B (gg 
v B (qiqj 



qq) 



<? B {qq -> 99) 



v B (.99 



<? B {qg 



99) 

qg) 



n .f 



ira Q 



ira s a s 



4 16m? 



27s 
+ 



27 



+ 



irat 



4 8m' 
9 I 27 ' 27s 



4m 



+ 



'8m 2 - 8mi\ 
~27s" + 27s 2 J 

4 8m 2 . 



Li 



9 9(m?s + mi) 



+ 



9s 



Li 



7ra 



^ 24 + 12s 



4m - 



m- 



+ 



+ r| log 



3s 3s 2 



+ Oil 



4m_ 



9 9(m~s + m_ 



+ 



4 
9 



5m; 
~9s~ 



1 + & 



Li 



8m? 



27(s + 2m 2 _ 



7TQ 



8 16mJ 

9 + 9s 



+ 



7ra s a s 
s 

2 



09 



8m. _ 



+ 



7ra 



32 



3s 
32mi 



+ 



8m? 8mi 
3s 



3s 2 



L 2 



27 ' 27(m?s + mi 



+ 



64m_ 



8m? 



7ra c 



S 
+ 



51m? N 
7 32m 2 



+ 



9.s 



2mf 



+ 



9 9m 2 
I s 

/ 8m 2 



27s 27(s - 2m 2 



log 



2m 2 -m 2 _ 



9s 



+ 



+ 



1 + /% 
8m 4 "■ 



9s 2 



+ 



2m?m?_ 



Li 



3 The explicit boundaries of these integrations are given in Appendix |b[ 
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with 



Li = Ioe 



Lq = lOE 



S + 2771,2 _ 

s + 2m 2 + s/3o 



S — 7771 — K 



S — 7771 + K 



Ami 



777_ = 777^ — 777^ 



a s = 5s/4vr 



Lo = lot 



La = lot 



S — 27T7 2 — 
S — 2777 2 + 



S + 777_ — K 
S + TO 2 , + K 



Ami 



/> — A / (s — 777? — 777?) 2 — 4777? 



9 m l 



Note that we have suppressed the theta- functions @(s — 4m 2 ) for the production thresh- 
olds. For identical particles in the final state [gluino-pair production or production of 
squarks with identical chirality and flavour] a statistical factor 1/2 has been taken into 
account. 

For the production of squark pairs or squark-gluino pairs only one initial state con- 
tributes at lowest order. For squark-antisquark and gluino pairs both gluon-gluon and 
quark-ant iquark initial states are possible. 

The total hadronic cross-sections are obtained by integrating the parton cross-sections 
in the usual way over the parton distributions fa in the proton/antiproton: 

<7{ij^Q,g) = J dx 1 dx 2 fi(x 1 )f j ( y X2)cr B ( y ij->q,g;s = x 1 x 2 S), (12) 

where the total centre-of-mass energy of the collider is denoted by \^S. 

To assess the relative weights of qq, qq, gg and qg final states in pp collisions at the 
Tevatron and in pp collisions at the LHC, the relative yields are shown for a typical set 
of mass parameters in Fig. ^ and Fig. [| The relative yields of squarks and gluinos in 
the final states depend strongly on the mass ratio mq/mg, for which we have chosen two 
representative values, 0.8 and 1.6. If squarks are lighter than gluinos, the valence partons 
give the dominant yield of squark-antisquark pairs/squark pairs at the Tevatron/LHC. 
By contrast, if the gluinos are the lighter of the two species, their production is the most 
copious. 



3 SUSY-QCD Corrections 
3.1 Virtual Corrections 

In this subsection we will present the virtual QCD corrections to the partonic production 
cross-sections of squarks and gluinos. The technical set-up of the calculation will be de- 
fined, including the renormalization of the ultraviolet (UV) divergences. The calculations 
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Figure 2: The relative yields of squarks and gluinos in the final states at the Tevatron. 



The mass ratio mq/m g is chosen to be 



a) 0.8 and (b) 1.6. Also shown are the leading 

0; 



parton contributions for (c) qq and (d) gg final states. Parton densities: GRV 94 



renormalization and factorization scale Q 



for squarks, Q 



for gluinos, and 



Q 



+ rrig)/2 for squark-gluino pairs. 



are carried out in the MS renormalization scheme, which requires a careful analysis of 
the Yukawa (qqg) coupling g s in higher orders. 

3.1.1 Technical Set-Up 

The calculation of the 0(a s ) corrections to the reactions (2)-(7) involves the evaluation 
of the virtual corrections, i.e. the interference of the Born matrix element [after ghost 
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Figure 3: The relative yields of squarks and gluinos in the final states at the LHC. The 
mass ratio mg/m g is chosen to be (a) 0.8 and (b) 1.6. Also shown are the leading parton 
contributions for (d) gg final states. Parton densities: GRV 94 | 20| ; renormalization and 

+ rrig) 1 2 for 



m, 



factorization scale Q = rriq for squarks, Q = m g for gluinos, and Q 
squark-gluino pairs. [Note that qq/qg final states can only be generated by qq/qg initial 
states so that diagram (c) is trivial and not shown.] 



subtraction] with the one-loop amplitudes. The corresponding differential cross-section is 
given by 



2^V 



K 



(t - p|)0 - vl) - vis 



H 2 s 



e([t- P 2 2 }[u-p 2 2 }- P 2 2 s) 



dtdu " lJ T(l-e) 

x e(s - 4m 2 ) 6(s + t + u - pi - pj) (M b M v * + M V M B *) . (13) 
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The virtual (one-loop) amplitudes include self-energy corrections, vertex corrections, and 
box diagrams. For the virtual particles inside loops we use the complete super symmetric 
QCD spectrum: gluons, gluinos, all quarks, and all squarks. In Fig. £| we present a set 
of typical virtual corrections. In (a) the gluino self-energy is given. The diagram with 
reversed fermion-number flow contributes due to the Majorana nature of the gluinos. The 
vertex corrections are exemplified by a gauge vertex (b) and a Yukawa vertex (c). Typical 
examples for the different types of box diagrams are depicted in (d), (e), and (f). 

The divergences in the virtual corrections are regularized by performing the calcula- 
tions in n = 4 — 2e dimensions. These divergences consist of ultraviolet^], infrared (IR), 
and collinear divergences [also called mass singularities] , and they show up as poles of the 
form e~ % (i = 1,2). Since the virtual amplitudes are contracted with the Born matrix 
elements, all loop momenta will be contracted with themselves or external momenta. This 
leads to a great simplification of the tensor integrals appearing in the one-loop corrections. 
These tensor integrals are evaluated in n dimensions by means of an adapted version of 
the standard Passarino-Veltman tensor integral reduction |[21|| , and they are expressed 
in terms of scalar integrals. The coefficients of these scalar integrals are finite, and they 
have to be calculated up to 0(e 2 ). The divergences are contained in the scalar integrals. 
We have calculated these, using two different techniques. One is based on the Feynman 
parametrization, the other proceeds via the computation of the absorptive part by apply- 
ing the Cutkosky cutting rules in n dimensions, followed by the use of dispersion-integral 
techniques to get the real part. All analytical manipulations were performed with the 
help of the symbolic computer program FORM ||22|| . 

The case of equal masses for squarks and gluinos is calculated separately, in view of 
additional singularities. Divergent terms of the form log[(m? — m?)/m?] lead to additional 
1/e poles. It has been checked explicitly that the final cross-sections are nevertheless 
continuous at this point of the mass-parameter space. 

The singularity structure of the scalar integrals in n dimensions can be summarized as 
follows. The non-vanishing scalar one-point and two-point functions only give rise to UV 
poles. The derivative of the on-shell two-point functions^, and the three- and four-point 
functions are UV finite and give rise only to IR and collinear singularities. IR poles appear 
when a massless particle is exchanged between two on-shell particles; collinear poles show 
up when a massless particle splits into two massless collinear particles. Double poles are 
generated only when IR and collinear singularities are present at the same time. 

The 75 Dirac matrix, entering the calculation through the quark-squark-gluino Yukawa 
couplings, is treated in the 'naive' scheme. In this scheme the 75-matrix anticommutes 
with the other 7 M -matrices. This is a legitimate procedure at the one-loop level for 
anomaly-free theories. 

We have excluded the top-squarks from the final states, as discussed in the intro- 
duction. However, to carry out the NLO calculation consistently, we have to take into 

6 As expected, no quadratic UV divergences are generated in softly broken supersymmetric models. 
7 In the top-stop loop contributing to the gluino self-energy we insert widths for top, stop, and gluino 
states. 
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Figure 4: A selected set of Feynman diagrams for the virtual corrections, (a) Gluino self- 
energy, (b) quark-quark-gluon vertex [gauge coupling], (c) quark-squark-gluino vertex 
[Yukawa coupling], (d) two-point boxes, (e) three-point boxes, and (f) four-point boxes. 



14 



account the top-squarks inside loops. For the sake of simplicity we take them to be mass 
degenerate with the other squarks. For the top quark we use the mass m t = 175 GeV. 
Thus, the final results will depend on two free parameters: the squark mass rriq and the 
gluino mass trig. 

The definitions of the invariant energies and momentum transfers, and the gauge 
choices for internal and external gluons are the same as in lowest order. Faddeev-Popov 
ghost contributions have therefore to be taken into account in the gluon self-energy and 
in the three-gluon vertex corrections. 



3.1.2 Renormalization of UV Divergences 

The UV divergences can be removed by renormalizing the coupling constants and the 
masses of the heavy particles. The external self-energies are multiplied by a factor 1/2 
to properly account for the transition from Green's functions to the S-matrix. For the 
renormalization of the QCD coupling constant one usually resorts to the MS scheme, 
which involves n-dimensional regularization, i.e. fields, phase space, and loop momenta 
are defined in n dimensions. The UV 1/e poles are subtracted, together with specific tran- 
scendental constants, at an arbitrary subtraction point Qr, the charge-renormalization 
scale. 

In supersymmetric theories, however, a complication occurs. In n ^ 4 dimensions the 



MS scheme introduces a mismatch between the number of gluon (n — 2) and gluino (2) 
degrees of freedom. Since this 0(s) mismatch will result in finite non-zero contributions, 
the MS scheme violates supersymmetry explicitly in higher orders. In particular, the 
Yukawa coupling g s , which by supersymmetry should coincide with the gauge coupling 
g s , deviates from it by a finite amount at the one-loop level. Requiring the physical 
amplitudes to preserve this supersymmetric relation, a shift between the bare Yukawa 
coupling and the bare gauge coupling must be introduced in the MS scheme: 



9s = 9s 



4tt V 3 2 



9s 



1 + 3, 



(14) 



which effectively subtracts the contributions of the false, non- supersymmetric degrees of 
freedom [also called e scalars]. 

The need for introducing a finite shift is best demonstrated for the effective [one-loop 
corrected] Yukawa coupling, which must be equal to the effective gauge coupling in an 
exact supersymmetric world with massless gluons/gluinos and equal- mass quarks/ squarks. 
For the sake of simplicity we define the effective couplings T eS (Q 2 ) and T cS (Q 2 ) in the 
limit of on-shell quarks/squarks and almost on-shell gluons/gluinos, with virtuality Q 2 <C 
to- = m 2 ; in this limit the couplings do not contain gauge-dependent terms. In the 
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MS scheme we find, after charge renormalization: 



MS : r eff (Q 2 ) 
f eff (Q 2 ) 



1 + 

4:71 



mi 



log 

2 & 



Q 2 



m- 



9s<l + 



4:71 



N 



log 



log 

2 s 



7 
6 

+ 



Cf 
2 



(15) 



(16) 



The singular term 1/e represents the combination 1/e — 7^ + log(47r). The remaining 
1/e poles are IR and collinear singularities. Inspecting T eS and T eS , it is easy to prove 
that the difference between the two effective couplings coincides with the shift in Eq. fll^ 



Taking into account this finite shift of the bare couplings in the MS scheme, both effective 
couplings become identical at the one-loop level. In this way supersymmetry is preserved 
and the MS renormalization becomes a consistent scheme. 



An alternative renormalization scheme is the modified Dimensional Reduction (DR ) 
scheme in which the fields are treated in four dimensions, but the phase space and loop 
momenta in n dimensions. In this scheme no mismatch between bosonic and fermionic 
degrees of freedom is apparently introduced and supersymmetry is preserved ab initio. At 
the level of the effective couplings T and T , this is reflected in the equalities 



DR : T eS (Q 



9s <1 



T ett (Q< 



— N 
4% 



1 , (l? 
- - log — 

e \m\ 



2 \ml 



(17) 
(18) 



As a result, both couplings are identical order by order. [It should be noted that the 
transition from the effective gauge coupling in MS to the gauge coupling in DR involves 
a well-known finite renormalization a s N/(247r) = a s /(87c) fl2T 



In the following we use the MS renormalization scheme, supplemented by the finite 
shift of the Yukawa coupling. In this way supersymmetry is preserved on the one hand, 
while on the other hand the definition of the strong gauge coupling corresponds to the 
usual Standard-Model measurements. 

Below we list the various renormalizations needed for the production of squarks and 
gluinos. In order to preserve the form of the Ward identity given in Eq. @, non-zero 
particle masses have to be renormalized in an on-shell scheme. We have opted for a real 
mass renormalization, involving the subtraction of the real part of the on-shell self-energies 
at the real-valued pole masses. In the case of squarks and gluinos, this is equivalent to 
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replacing the bare masses in the lowest-order propagators by 
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The parameters rriq and are the pole masses. 

As discussed above, the couplings are renormalized in the MS scheme, including the 
finite shift of the bare Yukawa coupling given by Eq. fll4|). This leads to the following 
replacements of the bare couplings in the LO expressions: 
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bare 
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The first coefficient /3 of the SUSY-QCD (3 function can be decomposed into a sum of 
contributions from light and heavy particles: 



0o 





2 




3 




+ 



2 2 1 

>-5-5< B ' + 1 ' 



/# + A?- 



(19) 



In addition to the poles, also some logarithms are subtracted in order to decouple the heavy 
particles [top quark, squarks, gluinos] from the running of a s (Q%). In this decoupling 
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scheme the Q\ evolution of the strong coupling is determined completely by the light- 
particle spectrum [gluons and rif — 5 massless quarks]: 



d\og(Q 



9s(Q 2 R )P(g s ) = -a 2 s (Q 2 R )P^. (20) 



R) 



The methods described above to renormalize the UV divergences lead to cross-sections 
that are UV-fmite. Nevertheless, there are still divergences left. The IR divergences will 
cancel against the contribution from soft-gluon radiation. The collinear singularities, 
finally, will be removed by applying mass factorization. These steps will be discussed in 
detail in the next subsections. 



3.2 Real-Gluon Radiation 

Two important aspects of the real-gluon radiation, the split-up of the phase space in soft- 
and hard-gluon regimes, and the renormalization of collinear divergences by means of 
mass factorization, will be discussed in detail in the following subsections. 

3.2.1 Matrix Elements and Phase Space 

In order to complete the NLO evaluation of squark and gluino production we need, in 
addition to the afore-mentioned virtual SUSY-QCD corrections, also the corrections from 
real-gluon radiation. They are obtained from the LO partonic reactions by adding a gluon 
to the final state: 



Qi+Qj - 


^qk+qi+g 


(21) 


9+9 - 


^qi+qi+9 


(22) 


Qi+Qj ~ 


^qi+qj+g 


(23) 


qi+qi - 


^9+9+9 


(24) 


9+9 - 


^9 +9+9 


(25) 


qi+9 - 


-> qi +9 +g- 


(26) 



Again, the momenta of the initial-state partons are denoted by ki,k 2 , while the particles 
in the final states carry momenta Pi,P2, and k 3 . The charge-conjugate final states, which 
are not given here explicitly, will be taken into account when the hadronic cross-sections 
are calculated. A representative set of Feynman diagrams, contributing to the real-gluon 
amplitude Ai R , is given in Fig. |5|. We display some diagrams for (a) squark-antisquark 
production in quark-ant iquark annihilation, (b) gluino-pair production in gluon fusion, 
and (c) squark-gluino production in quark-gluon collisions. The internal and external 
gluon lines, and the 75 Dirac matrix, are treated in the same way as before. 
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(c) 





Tnrmnr<- 




Figure 5: A representative set of Feynman diagrams corresponding to real-gluon radia- 
tion: squark-antisquark production (a), gluino-pair production (b), and squark-gluino 
production (c). 



To evaluate the squared real-gluon matrix elements [.M^l 2 we define the following 
kinematical invariants JT7| : 



s = (h + k 2 ) 2 
S3 = (h+p 2 ) 2 - ml 
t = (h-p 2 ) 2 
u = (kx -P2) 2 
ue = (h-Pi) 2 ~ m\ 



S5 
s 4 

t' 
u' 
u 7 



(P1+P2) 2 

(h+Pi) 2 ~ m 2 q 



(27) 



h) 2 
h) 2 



(h - Pi) 2 - m 2 d . 
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For the sake of convenience we will use the following additional invariants: 



S3g = S3 + rrt~ - m~ s ig = s 4 + m\ - m~ 
s Ax = s 4 + m\ -p\ = 2(k 3 p 1 ) 

t 1 =t-m 2 ~ t g = t-m 2 ~ 

u\ = u — m~ u g = u — m~ 

ue g = ue + m 2 g - m 2 ~ u 7g = u 7 + m\ - m§. 

The squared matrix elements must be evaluated in n = 4 — 2e dimensions up to 0(e 2 ). 

After performing the n-dimensional three-particle phase-space integration, we find for 
the double-differential distributions 



2 dV* S 2 £ ^ \{t- P l)(u-pl)-p 2 s 

S — A ' 



/i 2 S 



Q([t-p 2 2 ][u-p 2 2 }-p 2 2 s 



dtdu " u 2r(l-2 £ ) 

X 6(5 - 4m2) ( 2 S+^)i- £ @ ( 2 ^) / ^ E l-^l 2 ' ( 28 ) 

with 2(k 3 pi) = s + t + u—p 2 —p 2 > 0. The n-dimensional angular integration is given ex- 
plicitly by J dQ n = J* sin 1_2£ (^ 1 ) d6\ sm~ 2e (9 2 ) d9 2 [see Appendix To perform the 
angular integrations, we isolate the coefficients of the form (s') k (s") 1 . The variables s' and 
s" denote kinematical invariants of the list in Eq. (27), and k and I are integer numbers. 
One of those invariants should contain both integration variables (&i,9 2 ), whereas the 
second should depend only on 9\. This can be achieved by means of partial fractioning, 
exploiting the fact that only five of the kinematical invariants are in fact independent)^. 
The angular integrals are therefore of the form 

= / sin 1 - 2£ (9 1 )d9 1 / sin- 2£ (^ 2 )^ 2 (a + 6cos^i)- /£ (A + 5cos^i + Csin^ 1 cos^)^. 
Jo Jo (29) 

Explicit analytical expressions for these angular integrals can be found in Ref. |I7| . In 
Appendix |B| we demonstrate how the kinematical invariants can be expressed in terms of 
the angular variables. 



3.2.2 Soft- and Hard-Gluon Radiation 

To identify the IR singularities, the phase space for gluon radiation is split into two distinct 
regimes, one describing soft gluons and the other describing hard gluons. This separation 
can be defined by introducing a cut-off parameter A in the invariant mass s Ax = 2(k 3 pi), 
corresponding to the radiated gluon and one of the heavy particles in the final state. The 

8 The relevant relations between the various invariants are given in Appendix |b] for the process of 
squark-gluino production, which represents the most general case. 
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cut-off parameter is chosen so small that it can be neglected with respect to any other 
mass scale in the process. In terms of the single-differential distribution da/dt the split-up 
takes the form 

da R r™ , d 2 a R f A , d 2 a s , dPa H , , 

ds ^ 1717. = / ds ** 7777. + / ds ^ -7717.- ( 30 ) 



dt J Q dtdu J Q dtdu J A dtdu 

The first term on the right-hand side of Eq. (|30|) represents the regime of soft-gluon 
radiation. In this regime the momentum of the radiated gluon, ks, tends to zero and an 
eikonal approximation can be applied, i.e. neglecting k% whenever possible. In the limit 
&3 — > 0, the (2 — > 3) kinematics is reduced to (2 — > 2) kinematics, and the kinematical 
invariants of Eq. (27) take the form 

s 5 ^s 2(k 3Pl )^0 2 (hp 2 )^0 (31) 

u — > u G /u G g^ui/ug u 7 /u 7g — > tl/tg, 

while the remaining invariants are not affected. 

After integration over the angles and over the invariant mass s^ xi singular expressions 
of the form e~ l (i = 1,2) are generated. The double poles correspond to configurations 
where IR and collinear singularities coincide. When the single-differential soft-gluon dis- 
tribution do s jdt is added to the virtual corrections, the sum is IR-finite. This sum, 
however, is not free of divergences until the collinear singularities are removed by means 
of mass factorization. 

The second term on the right-hand side of Eq. ( pOf ) represents the regime of hard-gluon 
radiation. In this regime only collinear singularities occur, generated when the radiated 
gluon is collinear with one of the initial-state massless particles. They show up in |A^ fi | 2 
as terms proportional to 1/t' or 1/u', which behave as 1/[1 ± cos(^)] (i = 1,2) and lead 
to 1/e poles after the angular integration. Also these collinear singularities have to be 
removed by means of mass factorization. 

As a result of the split-up of the phase space, terms of the form log*(A/m 2 ) (i = 1, 2) 
occur in both the soft and hard cross-sections. They come from the same terms that 
generate the IR singularities. If soft and hard contributions are added up, however, any 
A dependence disappears from the cross-sections in the limit A — > 0. 

3.2.3 Mass Factorization 

The collinear divergences, generated by the radiation of gluons [or massless quarks] , have 
a universal structure. The partonic cross-sections cr^, which contain the collinear singu- 
larities, have the following form, factorized to all orders of perturbation theory: 

2 d 2 aij(s,t x ,u x , fi 2 ,e) f dxx f 1 dx 



dt x du x 



I — I — y^r Zi (a;i,Q|,/i 2 ,£:) r mi (x 2 ,Q|,/i 2 ,e) 
Jo xx Jo x 2 r— ' 



l,m 

.2 d &im{s> t XJ U x , Qp) 



x s 2 ~™^i"*^*^*'> (32) 
dt x du x 



21 



t x = -2(k 2 p 2 ) 



u. 



-2(kip 2 ) s = xix 2 s i x = x 2 t x u x = xxu 2 



The indices i-m characterize the initial-state partons. The universal splitting functions 
Tij, representing the probability of finding, inside the parent particle j at the scale Q 2 F) a 
particle i with fraction x of the longitudinal momentum, contain the collinear divergences. 
They can be absorbed into a redefinition of the parton densities at NLO [24|, in general 
called mass factorization. Since the subtraction point of the mass-factorization procedure 
is arbitrary, the splitting functions will depend on the factorization scale Qf- Adopting 
the MS mass-factorization scheme we can write to 0(a s ) 



Fij(x, Q 2 F ,fi 2 ,e) = SijS(l 



~ + log | f 



(33) 



with 1/s — 1/e — 7£ + log(47r) as before. The hard-scattering (reduced) cross-sections cr^ 
are free of collinear divergences. They depend, like the splitting functions, on the scale 
Q F . In NLO they have the form 
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(34) 
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The Altarelli-Parisi kernels Pij{x) [p5] are given by [Tf = 1/2] 
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(35) 



(36) 

(37) 
(38) 



The parameter 5 is related to the IR cut-off parameter A through relations of the form 
S = A/(s + t x ) or 5 = A/(s + u x ), as can be read off from Eq. fl3"3|) by solving the 5- 
functions in the LO distributions. For the mass factorization of the collinear divergences 
related to gluon radiation, only the diagonal splitting functions Yu are required. [In 
SUSY-QCD, additional splitting functions are realized in the final-state distributions at 
very high energies. For the sake of completeness they are collected in Appendix p.] 

After performing the mass factorization in this way, the final results for the virtual 
corrections plus gluon radiation are free of collinear divergences. 
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We will use the standard MS mass-factorization scheme in which most of the ex- 
perimentally determined parton densities have been parametrized. The transition to the 
DR scheme is non-trivial and involves a careful matching for gluon-initiated heavy-particle 
production. [This was first observed for top production within standard QCD [T7| and is 
not related to supersymmetry aspects.] 



3.3 Final States with an Additional Massless Quark 

In this subsection we shall discuss the partonic reactions that can only be realized in 
next-to-leading order of the SUSY-QCD perturbative expansion. These reactions involve 
final states with an additional massless (anti) quark. In such reactions, explicit particle 
poles show up inside the allowed phase space, requiring the careful isolation of on-shell 
squark and gluino production and a subsequent subtle subtraction procedure of the poles. 

3.3.1 Crossing 

The reactions that involve the radiation of a massless (anti)quark only contribute at NLO: 



g+qj — >q k +qi+qi (39) 

qi+g — >qk+qi+qj (40) 

qi+g — >qi+qj+qj (41) 

g+q t — >g+g+qi (42) 

qi+g — >g +g +g< (43) 

g+g — >qt +g +qi (44) 

qi+qj — ► q~k+g +q~i (45) 

Qi+Qj — ► Qk+g +qi- (46) 



Again, the momenta of the initial-state partons are denoted by k\,k2, while those of 
the particles in the final states are denoted by Pi,pz, and k 3 . In Fig. || we give a few 
selected Feynman diagrams for (a) the squark-antisquark-quark final state and (b) the 
gluino-gluino-quark final state. 

In fact, only the matrix element for process (45) requires a new calculation. The 
squared matrix elements for the other reactions can be related to subprocesses (21), (23), 
(24), (26), and a posteriori (45) by means of crossing. This involves the interchange of 
the particles with either the momenta k\ and k 3 (1 <->• 3) or k 2 and k 3 (2 <-> 3). 

The (1 3) crossing corresponds to the replacement ki <-> — ^3. In terms of the 
kinematical invariants of Eq. (27), this is equivalent to the exchanges 

s <-> t' s 3 /s 3g <-> ui/ug s 4 /s4g <-> u 7 /u 7g , (47) 
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whereas the other invariants are not affected. In view of the different sign of the quark 
momentum inside the spinor sum, the resulting squared matrix elements have to be mul- 
tiplied by a factor (—1). 

Analogously, the (2 <-> 3) crossing corresponds to the replacement k 2 <-» —k 3 . In terms 
of the kinematical invariants of Eq. (27), this is equivalent to 

s^u' sa/ssg^+ti/tg s 4 /s 4g <-> u 6 /u 6g , (48) 

whereas the other invariants are not affected. Again the resulting squared matrix elements 
have to be multiplied by a factor (—1). 

The double-differential distributions are defined by Eq. (p8[) . In parallel to real-gluon 
radiation, the integrand is cast into the appropriate form by means of partial fractioning 
before the angular integrals are performed. Since no IR divergences are generated, the 
split-up into soft and hard regimes is not needed. 

However, the distributions contain initial-state collinear singularities, which can be 
removed by means of mass factorization [Eq. (|34])1. The non-diagonal splitting functions 
are needed, together with the LO distributions with one gluon more or one gluon less in 
the initial state. At this point the difference in the degrees of freedom for the gluons and 
quarks in n dimensions starts to play a role, leading to finite contributions to the reduced 
cross-sections. 

3.3.2 On-Shell Intermediate Squark/Gluino States 

Quarks in the final state can be decay products of on-shell squarks (q — > gq) if rriq > m g , 
or of on-shell gluinos (g — * qq) if m g > rriq. Since these processes occur at lowest order 



p - - «• » V 




Figure 6: Selected set of Feynman diagrams for subprocesses that involve additional 
massless quarks in the final state. Squark-antisquark production (a), and gluino-gluino 
production (b). 
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and since the branching ratios for the decays are large, these channels are the dominant 
production channels for quarks in the final state. If the squarks and gluinos are off 
shell, cf. Fig. |, the production cross-section is suppressed by the strong coupling a s and 
classified as a higher-order correction. However, by inspecting Fig. || it is obvious that 
some of the higher-order amplitudes are smoothly connected with the Born amplitudes. 
Formally they are marked by singularities, l/(p 2 — m~) or 1/ (p 2 — m~) , if the momentum 
flow approaches the squark or gluino mass. These problems can easily be solved by 
introducing the non-zero widths of squarks/gluinos and regularizing in this way the higher- 
order amplitudes. After subtracting the Breit-Wigner pole contributions from the higher- 
order diagrams, which are already accounted for by the Born diagrams, the cross-section 
for the production of squarks and gluinos is defined properly and no double-counting 
occurs. 

To exemplify this procedure we restrict ourselves to the case in which the squarks are 
heavier than the gluinos so that the 'stable' final states, with respect to SUSY-QCD, are 
2-gluino final states. This example exhibits the full scope of subtleties inherent in the 
regularization and subtraction procedures. [The singularities generated in the other case 
in which the gluinos decay to squarks, are treated in a similar way] To be specific, we 
consider the subprocess qg — > ggq. 

The last diagram of Fig. |6|b gives rise to a particle pole if the squark momentum 
approaches the mq mass shell; the other diagrams correspond to continuum gg production. 
The pole is regularized by introducing the non-zero squark width Tg, substituting for the 
propagator the Breit-Wigner form 

1 1 

p 2 — m?~ p 2 — m 2 ~ + irrtqTg 

Denoting the on-shell resonance contribution to the matrix element, defined for p 2 = m~, 
by Aires, and the sum of the off-shell resonance contribution and the gluino continuum 
contribution by M. rem -, the squared matrix element can be decomposed into the resonance 
part and a remainder, 

\M\ 2 = \M res \ 2 + 2 Re [M* res M rem ] + \M rem \ 2 . (50) 

Integrating over the entire phase space of the ggq final state, the resonance contribution to 
the cross-section represents the qg Born cross-section [including the branching ratio T(q — > 
qg)/Tg], while the remainder is to be attributed to the 0(a s ) higher-order corrections to 
gg production!. 

The particle pole l/(q 2 — m~) associated with the other final-state gluino can be treated 
in a similar way. The presence of complex masses in the p and q propagators gives rise 
to real contributions from the interference of the two imaginary parts. This requires a 
careful treatment of the angular integrations [see Appendix [FJ]. Single poles in \Ai\ 2 of the 
form 1/ (p 2 — m ~) or 1/ (q 2 — m~) , corresponding to configurations with only one on-shell 

9 Technical details on the separation of the resonance contribution are deferred to the Appendices. 
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(49) 



propagator, give rise to principal-value integrals, resulting in finite contributions to the 
cross-sections. In these contributions we use a very small decay width for the squarks 
instead of the physical width Tg. The influence of the actual size of the squark width was 
found to be very small. 

Adding the 0(a s ) corrections to the process qg — > qg, the final result for the cross- 
section a(qg — > ggq) including all radiative corrections to 0(a s ), may therefore be written 
as: 

<t(Q9 -> 990) = °res(qg 990 [LO + NLO]) + Aa(qg -> ggq [rem.]), (51) 

with Aa denoting the interference term and the continuum contribution in Eq. flSOP . By 
definition we will attribute a res to qg production, but Aa to gg production. 

For squark-antisquark and squark-squark final states, similar procedures have to be 
followed if the gluinos are heavier than the squarks. For squark-gluino final states the 
subtraction procedure is required for nig > m g as well as for m g > nig. 

After all singularities have been removed, we end up with well-defined double-differen- 
tial distributions for the irreducible squark-antisquark, squark-squark, gluino-gluino, and 
squark-gluino final states. These irreducible final states include only those topologies in 
which the lightest of the coloured SUSY particles is not produced in on-shell decays of 
the heavier particle. 



4 Results 

4.1 Partonic Cross-Sections 

We first present the NLO SUSY-QCD results at the parton level for the production of 
squarks and gluinos in quark and gluon collisions. To classify the contributions it is con- 
venient to decompose the partonic cross-sections into scaling functions. In contrast with 
the double-differential cross-sections, these scaling functions for the total cross-sections 
can in general not be presented in analytic form. Nevertheless, for two kinematical limits, 
at high energies and close to the production threshold, compact analytical expressions 
can be derived. 



4.1.1 Scaling Functions 



The partonic cross- sect ions can be calculated from the double-differential distributions, 
discussed in the previous subsections, by integration over the Mandelstam variables t and 
S4. The exact boundaries for the integration can be found in Appendix [Bl For a detailed 
analysis of the partonic cross-sections we introduce scaling functions 
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The indices i, j = g,q,q indicate the partonic initial state of the reaction. As before, 
m = (\/pf + \/p|)/2 is the average mass of the produced particles. The centre-of-mass 
energy of the partonic reaction is absorbed in the quantity 77, which is better suited 
for analyzing the scaling functions in the various regions of interest. Note that we have 
identified the renormalization and factorization scales, Q = Qr = Qp, properly justified 
in the next subsection. For identical particles in the final state, i.e. gluino pairs or 
squark pairs with equal flavours and chiralities, we have taken into account the statistical 
factor 1/2. The scaling functions are divided into the Born term f B , the sum of virtual 
and soft-gluon corrections f v+s , the hard-gluon corrections f H , and the scale-dependent 
contributions /. In this context it should be noted that the log*(A/m 2 ) terms (i = 1,2) 
are removed from the soft-gluon corrections and added to the hard-gluon part. The hard- 
gluon corrections are therefore independent of the cut-off for A <C m 2 . 

The scaling functions for squark-antisquark production are displayed in Fig. [7|. Unless 
stated otherwise, we use rriq = 280 GeV, m g = 200 GeV, and m t = 175 GeV as mass- 
parameter input, representing an allowed mass configuration close to the present exclusion 
boundaries. In LO the only possible initial states are gluon-gluon (a) and quark-ant iquark 
states [with equal (b) and different (c) flavours]. The gluon-quark (d) channel is only 
realized at NLO. 

In Fig. |8| we present the scaling functions for squark-pair production. To lowest order, 
the final states are generated exclusively in quark-quark collisions [with equal (a) and 
different (b) flavours]. At NLO the gluon-quark (c) initial state starts to contribute. 

The scaling functions for gluino-pair production are displayed in Fig. |^. In LO this final 
state can only be produced in gluon-gluon (a) and quark-antiquark (b) reactions. Yet 
again, the gluon-quark (c) initial state is possible at NLO. The adopted mass configuration 
[with nig > rrig] allows an on-shell intermediate squark-gluino state, with subsequent 
decay of the squark into a gluino and a massless quark. As discussed in Section |3.3.2j , 
the associated singularity has to be subtracted in order to avoid double counting. After 
the subtraction is performed, a remaining, though integrable, singularity shows up in 
the scaling function fg at the threshold for squark-gluino production. This integrable 
singularity is regularized by using a small non-zero squark width. 

The scaling functions corresponding to gluino-squark production are given in Fig. |H]. 
Solely the gluon-quark (a) initial state contributes at LO. All other initial states, i.e. 
gluon-gluon (b), quark-antiquark (c), and quark-quark (d), appear only at NLO. The 
singularities associated with on-shell squark-(anti) squark intermediate states are sub- 
tracted, leaving behind an integrable remnant. 

The comparison of the various scaling functions reveals that contributions involving at 
least one gluon in the initial state and at least one gluino in the final state are dominant. 
This is a straightforward consequence of the large colour charge of particles in the adjoint 
representation. 
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Figure 7: The scaling functions [Eq. (|52|)1 for squark-antisquark production from (a) gg, 
(b) qq, (c) q'q, and (d) gq initial states. Flavours and chiralities of the squarks are summed 
over. Mass parameters: rriq = 280 GeV, m- g = 200 GeV, and rrit = 175 GeV. 
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Figure 8: The scaling functions [Eq. (0)] for squark-pair production from (a) qq, (b) q'q, 
and (c) gq initial states. Flavours and chiralities of the squarks are summed over. Mass 
parameters: rriq = 280 GeV, m~ g = 200 GeV, and mt — 175 GeV. 
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Figure 9: The scaling functions [Eq. ([52|)] for gluino-pair production from (a) gg, (b) 
qq, and (c) gq initial states. The (integrable) singularity in (c) is the result of on-shell 
intermediate squark-gluino states; this singularity is regularized by using a small non-zero 
squark width. Mass parameters: nig = 280 GeV, m g = 200 GeV, and m t = 175 GeV. 
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Figure 10: The scaling functions [Eq. (|52])] for squark-gluino production from (a) gq, (b) 
gg, (c) qq, and (d) qq initial states. Flavours and chiralities of the squarks are summed 
over. The (integrable) singularities in (b) and (c) are caused by on-shell intermediate 
squark-antisquark states, in (d) by on-shell intermediate squark-squark states; these sin- 
gularities are regularized by using a small non-zero squark width. Mass parameters: 
nig = 280 GeV, m- g = 200 GeV, and m t = 175 GeV. 
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Figure 11: The variation of the scaling functions [Eq. (0)] for qq — > gg over a squark-mass 
interval around rriq = mf rrig = 200 GeV and (a) rriq = 175 GeV, (b) rriq = 200 GeV, (c) 



nig = 225 GeV. 



Noteworthy is the squark-mass dependence of the cross-section for gluino-pair pro- 
duction from quark-antiquark annihilation. This dependence is exemplified in Fig. |TT 
for the scaling functions, using the mass parameters m g = 200 GeV, as before, and 
rrig = 175, 200, 225 GeV. In the region of almost mass-degenerate squarks and gluinos, 
the maximum of the LO scaling function f B decreases with increasing squark mass. By 
contrast, the virtual and soft corrections f v+s increase for small and intermediate ener- 
gies [rj < 10], as is also evident from Fig. |9|b. The hard-gluon corrections f H are nearly 
independent. The ratio of the higher-order correction over the lowest-order cross-section 
will therefore vary rapidly for gluino-pair production in the range where gluino and squark 
masses are of the same order. 



4.1.2 Threshold Region 

The energy region near the production threshold is the base for an important part of the 
contributions to the hadronic cross-sections. This region is characterized by the small 
velocity f3 of the produced heavy particles in their centre-of-mass system ^ < 1]. Two 
sources of large corrections can be identified in this threshold region, the leading terms of 
which can be calculated analytically at NLO. These analytical expressions provide pow- 
erful checks of the numerically integrated NLO corrections for arbitrary mass parameters. 

First of all, the exchange of (long-range) Coulomb gluons between the slowly moving 
massive particles in the final state [see Fig. [12|a] leads to a singular correction factor 



~ 7ra s //9, which compensates the LO phase-space suppression factor (3 [p^j . The scaling 
function f v+s therefore tends to a non-zero constant at threshold. It should be noted, 
however, that the screening due to the non-zero lifetimes of the squarks/gluinos reduce 
this effect considerably. 
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Figure 12: Generic diagrams leading to (a) the Coulomb singularity, (b) the large thresh- 
old logarithms, and (c) the high-energy plateau. The solid external lines represent glu- 
ons/(anti) quarks, the dashed ones respresent gluinos/(anti)squarks. 



Secondly, as a result of the strong energy dependence of the cross-sections near thresh- 
old, large positive "soft" corrections ~ log J (/3 2 ) (i = 1,2) are observed in the initial-state 
gluon-radiation contribution [see Fig. H|b]. They can in fact be resummed [^7, 2q| . In 
our analysis we will, however, stick to the strict NLO corrections. 

Near threshold the scaling functions can be expanded in (3, leading to the following 
analytical expressions [suppressing 0(/3)]: 
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For squark-antisquark final states the scaling functions for equal and different flavours 
are identical (f q 'g = f q q) near threshold. 

Squark— Squark: 
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Gluino— Gluino: 
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The LO and NLO threshold cross-sections for gluino-pair production from quark-ant iquark 
annihilation vanish if the squarks and gluinos are mass degenerate. This follows from the 
destructive interference between the three LO diagrams. In Fig. [11] this phenomenon is 
clearly visible. 
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Note that for squark-gluino production the leading corrections factorize only partly in 
terms of the LO cross-section. This in contrast to the other production processes, which 
involve equal-mass final-state particles. 



4.1.3 High-Energy Region 

At high energies the NLO partonic cross-sections can asymptotically approach a non-zero 
constant, rather than scaling with 1/s as the LO cross-sections: 

a ~ ^ (LO) (57) 

s 

a ~ f?L (NLO). (58) 
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This is caused by almost on-shell, soft gluons in space-like propagators, associated with 



hard-gluon/quark radiation [see Fig. 0c]. Since the splitting probabilities q — > qg and 
g — ► gg are scale-invariant mod. logarithms, the size of the NLO cross-section is set by 
the centre-of-mass energy of the subprocess induced by the virtual gluon (marked "Born" 
in Fig. |l2|c). For the dominant contributions this centre-of-mass energy is of the order 
of the squark/gluino masses, i.e. not far above the threshold. It should be noted that 
these high-energy plateaus only have a marginal influence on the hadronic cross-sections, 
as the main part of the contributions originates from the partonic energy region near the 
production threshold. 

Exploiting the factorization in the transverse gluon momentum at high energies |29| , 
the high-energy scaling functions can be determined analytically. 



Squark— Antisquark: 



fH = ^9_ fH= _21W_ m 
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The ratio of the f gg and f gq scaling functions is given by 2N : Cp = 9 : 2. This ratio 
corresponds to the probability of emitting a soft gluon from a gluon [~ N, 2 sources] or 
a quark (~ Cp, 1 source). 

Squark— Squark: The squark-pair production cross-section does not exhibit a high- 
energy plateau, since no space-like gluon-exchange diagrams are possible at NLO. 

Gluino— Gluino: 

fH _ 1949 rH_ 1949 

fgg ~ 160tt Uq ~ 240tt- 
The ratio of the f gg and f gq scaling functions is given again by 2N : Cp = 9 : 2. 

Squark— Gluino: 

fH _ 517 H = rH_ = 517 fH _fH _ 517 , . 

iqg 864tt Jqq iq " q 1944tt Jqq Jq ' q 972tt 1 ; 
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<Q9- ^ JQQ -JQ'Q- ^ JQQ -JQ'Q- g^- 

These (simple) high-energy limits are derived for equal squark and gluino masses; the 
results for arbitrary masses are very complicated and are therefore not given explicitly. 
The ratio of the f qg , f qg , and f qq scaling functions is given by N : Cp : 2Cp = 9:4:8, 
irrespective of the precise values for the squark and gluino masses. 
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4.2 Hadronic Cross-Sections 



Finally we discuss in this subsection the hadronic cross- sect ions for the production of 
squarks and gluinos. The analyses are performed for the Fermilab pp collider Tevatron 
with a centre-of-mass energy of VS = 1.8 TeV, and for the CERN pp collider LHC with 
a centre-of-mass energy of vS = 14 TeV. In analogy to the experimental analyses, we 
consider the following four hadronic production processes 

pp/pp — > qq (62) 

pp/pp — > qq,qq (63) 

pp/pp -> gg (64) 

pp/pp _> g#,g<?. (65) 

As before the chiralities and flavours of the squarks (e.g. UL,dn) are implicitly summed 
over. Yet stop production is not taken into account; these final states will be analyzed in 
a forthcoming report ||. The charge-conjugate final states are now properly taken into 
account in the reactions (^) and flBSj). From now on we will refer to these two hadronic 
reactions, for simplicity, as qq (squark-squark) and qg (squark-gluino) production, re- 
spectively. 

Various aspects of the above reactions will be presented in detail. First of all, the 
dependence of the cross-sections on the renormalization/factorization scale and parton 
densities is investigated. Then the NLO corrections are studied for a default choice of 
scale and parton densities. We discuss the NLO effects on the differential distributions 
with respect to the rapidity y and transverse momentum p t of one of the outgoing particles. 
Finally we describe the NLO effects on the total cross-sections and their implications on 
the experimental search for squarks and gluinos. 



4.2.1 Scale and Parton-Density Dependence 

The total hadronic cross-sections are obtained by convoluting the partonic cross-sections 
with the relevant parton densities [in the proton or antiproton] : 



a(S,Q 2 ) = V f dx t I dx 2 f^(x 1 ,Q 2 )f^(x 2 ,Q 2 )a ij (x 1 x 2 S,Q 2 ) 



>2 - 



The partons i and j carry fractions x\ and x 2 of the original momenta of the hadrons h\ 
(= p) and h 2 (= p/p), respectively. The integrations in Eq. ( p^| ) are performed numerically, 
using the Monte Carlo integration routine VEGAS [0 . The parton densities have been 
extracted from a large variety of experiments and are available in various parametrizations. 
To estimate the associated uncertainty, we compare the results for a sample of three 
different parametrizations. 
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As a first step, we present for reactions (p^j -QBq) the dependence of the total cross- 



section on the renormalization and factorization scale Q = Qr = Qp in Figs. 13 and 



14| . We have checked that the separate variation of the factorization scale Qp and the 
renormalization scale Qr leads to a variation of the next-to-leading order cross-section 
(roughly) within the band generated by Q. We can therefore keep the discussion to this 
simplified case without loss of generality. The scale is restricted to Q < 1 TeV, as the 
parton densities are not available beyond this value. 

The results for the Tevatron are given in Fig. using the mass parameters nig = 
280 GeV, m g = 200 GeV, and m t = 175 GeV. For a consistent comparison of LO and 
NLO results, we calculate all quantities [a s (Q 2 R ), the parton densities, and the partonic 
cross-sections] in LO and NLO, respectively. In LO we use the parton densities of GRV94 
and CTEQ3 |HJ with the corresponding QCD couplings. At NLO this is compared 



with the parton densities of GRV94 0, CTEQ3 JEJ, and MRS(A') In LO the scale 



dependence is steep and monotonic. Changing the scale from Q = 2m to Q = m/2, the 
cross-section increases by 100-120%. In NLO the scale dependence is strongly reduced, to 
about 40-50% in this interval. At the same time the cross-section is significantly enhanced 
at the central scale (Q = m). The uncertainty originating from the various parametriza- 
tions of the parton densities in NLO amounts to < 10% at the central scale. An exception 
is the squark-squark cross-section, which is dependent on the badly determined sea-quark 
distribution. 



The results for the LHC are given in Fig. |TJ, using the mass parameters trig = 600 GeV, 
trig = 500 GeV, and m t = 175 GeV. Here, too, a strong reduction of the scale dependence 
is observed by taking into account the NLO corrections, as well as a clear enhancement of 
the cross-section at the central scale. In the interval between Q = m and Q = m/2 the LO 
cross-section increases by about 35%, whereas the variation for the NLO results is reduced 
to only 5-10%. At the LHC the dependence on the factorization scale is very weak and 
the residual scale dependence is dominated by a s [in contrast with the Tevatron where 
fairly large x values in the parton distributions give rise to a stronger factorization-scale 
dependence] . The uncertainty due to different parametrizations of the parton densities in 
NLO amounts to < 13% at the central scale. This is a result of the prominent role played 
by the gluon densities at the LHC. 

In conclusion, for all four reactions and for both hadron colliders the scale dependence 
is reduced by a factor of 2.5-4 when the theoretical predictions are improved by taking 
into account the next-to-leading order SUSY-QCD corrections. Even a broad and shallow 
maximum develops at scales near one third of the central scale. In the following we will 
adopt GRV94 parton densities and take Q = m as the default scale; this results in a 
conservative estimate of the cross-sections. 
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For charm and bottom quarks we use the earlier GRV distributions |31| 
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10 _1 1 Q/m s 10 1 1 2Q/(m~+ mi ) 



Figure 13: The dependence on the renormalization/factorization scale Q of the LO and 
NLO cross-sections for (a) squark-antisquark, (b) squark-squark, (c) gluino-gluino, and 
(d) squark-gluino production at the Tevatron {\fS = 1.8 TeV). Parton densities: GRV94 
(solid), CTEQ3 (dashed), and MRS (A') (dotted). Mass parameters: m~ q = 280 GeV, 
m-g = 200 GeV, and m t = 175 GeV. 
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Figure 14: The dependence on the renormalization/factorization scale Q of the LO and 
NLO cross-sections for (a) squark-antisquark, (b) squark-squark, (c) gluino-gluino, and 
(d) squark-gluino production at the LHC (\/S = 14 TeV). Parton densities: GRV94 
(solid), CTEQ3 (dashed), and MRS(A') (dotted). Mass parameters: m~ q = 600 GeV, 
m- = 500 GeV, and m t = 175 GeV. 
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4.2.2 Differential Distributions: Transverse Momentum and Rapidity 



The hadronic differential distributions are obtained by convoluting the partonic double- 
differential distributions with the relevant parton densities. We consider the distribu- 
tions with respect to the transverse momentum p t and rapidity y of one of the outgoing 
squarks/gluinos. The corresponding double-differential hadronic distribution is given byf 1 ] 



111(7 2 Pt S [ dx 2 xJ^(x u Q 2 )x 2 f^(x 2 ,Q^ !!2fT >'' r]r2S ' Q 



2^ 



dptdy . ' _J X - J x - 1 J dtdu 



(67) 



The definitions of pt and y and of the integration boundaries can be found in Appendix |B|. 
Squarks and antisquarks are not distinguished in the final state. So, for qq, qq, and gg 
final states we have to add the distributions with respect to both final-state particles. 



The rapidity distributions, as presented in Figs. 16 and 18, are defined as the sum of the 



contributions of positive and negative rapidity. The distributions are normalized to unity. 



In Fig. [15| the normalized p t distribution is given for the Tevatron. The input mass 
parameters are rriq = 280 GeV, trig = 200 GeV, and m t = 175 GeV. For the scale Q and 
for the parton densities we take the default settings [Q = m and GRV94]. In order to 
perform a consistent comparison of LO and NLO results, the parton densities are used 



with the associated values for a s . In Fig. [16] the normalized y distribution is shown. The 



corresponding distributions for the LHC can be found in Figs. [17] and |18[ elaborated for 
rrig = 600 GeV, m g = 500 GeV, and m t = 175 GeV. 

The normalized p t distributions are hardly affected by the transition from LO to NLO. 
In total, the NLO corrections render the distributions a little softer. This is caused by 
the fact that for high-energetic massive particles the probability of losing energy through 
radiation is large. The normalized rapidity spectra in LO and NLO are identical for all 
practical purposes. 

In conclusion, the properly normalized distributions of the squarks and gluinos in 
transverse momentum and rapidity are described quite well by the lowest-order approxi- 
mation. 



11 The preferred scales Q are the transverse masses; however, for convenience we have chosen the particle 
masses, as for the total cross-sections. 
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Figure 15: Normalized transverse-momentum distributions in LO (dotted) and NLO 
(solid) at the Tevatron (\^S =1.8 TeV). Parton densities: GRV94, with scale Q = m; 
mass parameters: rrig = 280 GeV, rrig = 200 GeV, and m t — 175 GeV. 
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Figure 16: Normalized rapidity distributions in LO (dotted) and NLO (solid) at the 
Tevatron (\^S = 1.8 TeV). Parton densities: GRV94, with scale Q = m; mass parameters: 
rrig = 280 GeV, m- g = 200 GeV, and m t = 175 GeV. 
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Figure 17: Normalized transverse-momentum distributions in LO (dotted) and NLO 
(solid) at the LHC (v 7 ^ = 14 TeV). Parton densities: GRV94, with scale Q = m; mass 
parameters: nig = 600 GeV, nig = 500 GeV, and m t — 175 GeV. 
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Figure 18: Normalized rapidity distributions in LO (dotted) and NLO (solid) at the 
LHC (\^S = 14 TeV). Parton densities: GRV94, with scale Q = m; mass parameters: 
nig = 600 GeV, m g = 500 GeV, and m t = 175 GeV. 
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4.2.3 Total Cross-Sections for Squark and Gluino Production 



fC-factors: To facilitate the quantitative comparison of LO and NLO cross-sections we 
define the ratio 

K = ctnlo/vlo, (68) 

usually referred to as the if-factor. For consistency, the cross-section a^o (&nlo) is 
calculated for all entries taken at leading (next-to-leading) order, i.e. couplings, parton 
densities, and parton cross-sections. 

In Fig. [19] we present the .fT-factors at the Tevatron for the reactions fl62|)-(|65|). For 
the scale Q and the parton densities we take the default settings [Q = m and GRV94]. 
In Fig. |19| a the i^-factors are displayed for a gluino mass m g = 200 GeV and for squark 
masses in the range 150-400 GeV, whereas in Fig. [l9| b the role of the squarks and gluinos 
is interchanged. The corrections strongly depend on the process. With the exception 
of the squark-pair production process, which is rather unimportant at the Tevatron, all 
processes are subject to large, positive corrections between +10% and +90%. The K- 
factors for squark final states are almost mass-independent. By contrast, the i^-factors 
for the (dominant) final states that involve at least one gluino exhibit a strong mass 
dependence. The particularly strong mass dependence of the gluino-pair .fT-factors for 
almost mass-degenerate squarks and gluinos is a direct consequence of the phenomena 
described in Section |4.1.1| for the scaling functions, since a large part of the hadronic 



cross-section originates from the quark-ant iquark channel. For a fixed gluino mass and 
increasing squark masses the LO cross-section decreases, whereas the virtual corrections 
increase. This leads to the observed increase of the i^-factor. If the squark mass is kept 
fixed and the gluino mass is increased, the reverse is observed. 



In Fig. |20| the if-factors are presented for the LHC. The input is the same as before, 
except for the fact that we consider fixed ratios of the squark and gluino masses: irtq/mg = 
0.8, 1.2, 1.6 and 2 [Figs. ^ a - d ]. Again, the corrections are positive and in general large, 
between +5% and +90%. The mass dependence and the absolute size of the .fT-factors 
for squark final states are moderate. By contrast, for final states involving gluinos the 
corrections are substantially larger and exhibit a strong mass dependence. The influence of 
the squark mass on the gluino-pair cross-section is less pronounced than for the Tevatron, 
since the gluon-gluon initial state yields the dominant contributions. 



Total cross-sections: The absolute size of the total hadronic cross-sections plays a 
crucial role in the experimental analyses. As long as no squarks and gluinos are discovered, 
the exclusion limits of the masses are derived by comparing the experimental data with 
the expected rates based on the theoretical predictions of the cross-sections. If squarks 
and gluinos are discovered, the masses of the particles will be determined experimentally 
by the same method. 

In Fig. [H] the total cross-sections are given for the Tevatron. The NLO results are 
based on the default settings: renormalization/factorization scale Q = m and GRV94 par- 
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Figure 19: The f^-factors [Eq. (p8|)1 for the Tevatron [\[S = 1.8 TeV). Parton densities: 
GRV94, with scale Q = m; mass parameters: m t = 175 GeV and (a) m g = 200 GeV, (b) 
nig = 200 GeV. 
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Figure 20: The iT-factors [Eq. (H)] for the LHC (VS = 14 TeV). Parton densities: 
GRV94, with scale Q = m; top-quark mass: m t = 175 GeV. 
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ton densities. These cross-sections are compared with the LO parametrizations adopted 
in the experimental analyses [11], [12| until recently [EHLQ parton densities, with Q equal 



to the partonic centre-of-mass energy^. Over the full mass range covered by the Teva- 
tron, the net effect of the NLO corrections is to raise the derived lower mass bounds by 
+ 10 GeV to +30 GeV. 



In Fig. the total cross-sections are given for the LHC, using the default settings 
and a representative range of squark and gluino masses. Over the full mass range covered 
by the LHC, the cross-sections are increased by the NLO corrections, leading to a shift 
in the associated particle masses in the range between +10 GeV and +50 GeV. 



4.2.4 Implications for Experimental Searches 

The precise knowledge of the cross-sections at next-to-leading order SUSY-QCD has a 
profound impact on the experimental analyses: 

(i) The renormalization/factorization scale dependence is reduced by roughly a factor 
of 2.5-4 compared with leading-order calculations, and the theoretical predictions of the 
cross-sections are stable. Taking for the renormalization/factorization scale Q the average 
mass m of the produced particles results in a conservative estimate for the cross-sections 
at NLO. 

(ii) The NLO corrections are large and positive at the central scale Q = m. The NLO 
corrections must therefore be included in the analyses to obtain adequate theoretical 
predictions for the total cross-sections, as required for deriving experimental mass bounds 
or measuring the squark and gluino masses. 

(in) The shape of the differential distributions in transverse momentum and rapidity 
of one of the outgoing squarks or gluinos is hardly affected by the NLO corrections. 

(iv) The NLO cross-sections raise the present lower mass bounds for squarks and 
gluinos derived from Tevatron data by +10 GeV to +30 GeV. 



12 The scale choice Q = s/s is only legitimate in LO; in NLO a hadronic scale, e.g. a particle mass, is 
required by the renormalization group. 
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Figure 21: The total cross-section for the Tevatron (y/S = 1.8 TeV). NLO (solid): GRV94 
parton densities, with scale Q = m; compared with LO (dashed): EHLQ parton densities, 
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5 Conclusions and Outlook 



In this report we have presented the next-to-leading order SUSY-QCD corrections for 
the production cross-sections of squarks and gluinos at the hadron colliders Tevatron and 
LHC. By reducing the scale dependence of the cross-sections considerably, the quality 
of the theoretical predictions is substantially improved compared with the lowest-order 
calculations. The NLO cross-sections provide a solid basis for experimental analyses of 
squark and gluino mass bounds/measurements at hadron colliders. 

So far the calculations have been performed for mass-degenerate squarks associated 
with the five light quark flavours. Generally, small mass differences between the L and R 
squark states for a given flavour are suggested by supergravity-inspired parametrizations 
of low-energy supersymmetry. Due to the large top-Higgs Yukawa coupling, however, the 
assumption of mass degeneracy is expected to be strongly broken for the two stop states t\ 
and t 2 , mixtures of the L and R chirality states II and Ir. In these cases the NLO SUSY- 
QCD cross-sections require the extension of the theoretical analysis to different left- and 
right-handed couplings of the quarks to squarks and gluinos, adding to the complexity 
of the calculation in a non-trivial way. For final-state stop particles this analysis is in 
progress and will be completed in due time. 



Note 

The Fortran codes of the NLO cross-sections can be obtained from hoepker @ x4u2.desy.de, 
spira @ cern.ch, or wimb @ lorentz.leidenuniv.nl. 



Acknowledgements 

We have benefited from discussions with K.I. Hikasa, M. Kramer, and W.L. van Neerven. 
Special thanks go to our experimental colleagues S. Lammel and M. Paterno for valuable 
comments on the Tevatron data and useful suggestions for experimentally convenient 
parametrizations of the theoretical results worked out in this report. 



48 



A Fermion-Flow Diagrams in SUSY-QCD 



The field-theoretic components of supersymmetric QCD are quarks/squarks and glu- 
ons/gluinos. The Majorana character of gluinos renders the evaluation of Feynman dia- 
grams involving these particles somewhat cumbersome. However, a simple prescription 
has recently been proposed in Ref. [34|, which allows an easy and fail-safe evaluation of 
the diagrams. It involves the definition of a continuous fermion-flow line, which in general 
does not coincide with the flow of the fermion number [given by the direction of the Dirac 
propagator line]. The amplitude must be evaluated along the fermion-flow according to 
the standard rules. In this way also fermion-number-violating processes like qq — > qq can 
be described in a straightforward fashion. 

This method is equivalent to the usual evaluation of SUSY-QCD diagrams if the 
analytic form of the propagators and vertices involving fermions is adjusted properly. The 
analytic expressions associated with these propagators and vertices are given explicitly in 
Fig. Indices of the fundamental /adjoint representation of colour SU (3) are denoted by 
i,j/a-c, the generators of the fundamental representation by t a = X a /2 (a = 1, 8), and 
the structure constants by f abc . The index \x is a Lorentz index, a and (3 are Dirac indices. 
The coupling constants g s and g s are the gauge and Yukawa couplings, respectively. The 
fermion-flow is defined in the diagrams by an additional directed line, parallel to the 
propagators. All other SUSY-QCD vertices and propagators are not affected by defining 
the fermion-flow lines. The flavour of the quarks and squarks are conserved in the SUSY- 
QCD interactions. 



B Kinematics and Phase Space 

The kinematics and phase space for the production processes of a pair of squarks or a pair 
of gluinos are the same as for top-antitop production |F7| . However, for the production of 



squark-gluino final states, the masses of the heavy particles are different. The kinematical 
and phase-space relations must therefore be derived for this more general case. Identifying 
squark and gluinos masses, we recover the formulae for squark and gluino pairs. 



B.l Partonic Processes 

We shall discuss the partonic process 

q(h) + g(h) — q( Pl ) + ~g(p 2 ) [+g(k 3 )} (69) 

for illustration. All particles are on shell, i.e. k\ = k\ = 0, p\ = m~, p\ = m| and 
k\ = 0. The kinematics of the radiative process is characterized by the invariant variables 
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Figure 23: Analytic expressions for the SUSY-QCD propagators and vertices that involve 
a specific fermion-flow, defined by the line parallel to the propagators. 

introduced in Eq. (27); they are related by energy-momentum conservation: 



5 4 = S + tg + Ml 

5 5 = S + t' + «' 

u 7 = —s — u„ — u'. 



S 3 = S + U 6 + U 7 
Uq = —S — t g — t' 



The total partonic cross-section is obtained by integrating the double-differential cross- 
section 



r4 r s 
a = / dt g 

Jt~ Jo 



dsA 



d 2 a 

dt„ ds^ 



(70) 
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within the limits 



s + m| — to- 



1 



to- — to?) 2 — 4m?m? 



s + £ + m s - to^ + 



TO?S 



The invariant energy of the qg subsystem in the final state is characterized by the variable 
s 4 , the square of the momentum transfer from the gluon (quark) in the initial state to the 
gluino in the final state by t = t g + m? (u = u g + to?). 



The double-differential cross-section in Eq. (|70|) can be derived by integrating the 
general four-fold differential cross-section over the angles 9^ defined in the centre-of-mass 
frame of the qg subsystem. In this particular frame the n-dimensional momenta are given 
by 

(71) 



h = 


(i0i,O, . 


., 0, 0, \p\ sinijj, \p\ cos ^ — u> 2 ) 


k 2 = 


(102,0, . 


.,0,0,0,ii; 2 ) 


h = 


(l^3,0, . 


., 0, W3 sin #1 sin 6*2, 103 sin 9\ cos #2, 103 cos^i) 


Pi = 


(#1,0,. 


. . , , — w 3 sin 9\ sin 6*2 , —103 sin 6*i cos 9 2 , —103 cos 9i 


P2 = 


(£2,0,. 


.., 0, 0, p sin^, \p\ cos^). 



The energies lOj, E^ the momentum and the angle ip are related to the ^-independent 



invariants s, t g , u g , and s 4 introduced earlier: 
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t 9 s 4s - s(w s + 2m? 
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4m? s 



Using these relations, all the invariant variables defined in Eq. (27) can be expressed in 
terms of s,t g ,u g , s 4 ,#j. For certain combinations of invariants, it is convenient to define 
the 0-axis with respect to k\ or p 2 (instead of k 2 ). Writing the n-dimensional angular part 
of the phase-space element as dVL n = sin 1_2£ (6 l i) d9\ sin _2e (#2) d9 2 , the double-differential 
cross-section is obtained from the matrix element A4 R in the following way: 



2„R 



d 2 a 



dt g ds A 
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&(tgU g - m 2 g s) Q(s - [ffiq + m g ] 2 ) 
2,— ®( s *) J dQ n Y,\M R \ 2 . (73) 
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Figure 24: Examples for on-shell intermediate resonance states: (a) gluino-pair interme- 
diate state for squark-gluino final states; (b) squark-antisquark intermediate state for 
squark-gluino final states; and (c) squark-gluino intermediate state for gluino-pair final 
states. 



In the limit rriq = rrig the well-known expressions of Ref. [17 are recovered for the equal- 
mass case. The 9-functions represent the requirements of positive energies and cos 2 ip < 1. 
They translate into the integration boundaries of the s 4 and t g integrals in Eq. (|70|) . 

As described in Section |3.3.2j , on-shell resonance production of intermediate particles is 
possible. For squark-gluino production, this must be analyzed carefully in the kinematical 
range where the three-particle final state gqq approaches the two-particle final state gg 
(if rrig > rriq) or qq (if rriq > rrig). This gives rise to singularities in l/s 3 or l/s 49 , as is 
evident from Fig. if the momentum flow approaches the q or g mass shells. 

After exchanging the order of the integrations, the total cross-section corresponding 
to the diagram in Fig. E3a can be written as 



'4 



a = dsi I dt g = / ds 4 — ^- (74) 

'0 Jtg(si) dtgdsn J s 4g 



i?-m?-20^f tf( Si ) = - S -^- ± 2 ^/(s - s ig ) 2 - 4sm§. 



The singularity in s<± g for m g > rrig can be regularized by inserting the non-zero gluino 
width and introducing the Breit-Wigner form, i.e. + wi|r~], for the (absolute) 

square of the propagator. 

Inserting the identity /(s 4£ ,) = [/(s 43 ) — /(0)] + /(0), the part of the cross-section 
related to /(0) corresponds to the (LO) production of an on-shell intermediate gluino and 
subsequent (LO) decay of this gluino into a squark and antiquark, since 



/(0) = 



B_ m g^g ^-g^q 

~ 9 ' 9 IT Tg 



and rrigTg/[s\ g + m?r?] — > 7r5(s 43 ) for small T g . This part is already accounted for by the 
lowest-order gg cross-section and has to be subtracted from the gqq cross-section. After 
the formal expansion of /(s 49 ) around s 49 = 0, the remaining part 



Act = / ds A (75) 
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is finite and well denned as a principal-value integral, since Sig/[s\ g + TTifr-] — > V(l/s4 g ) 
for small T g . 

For TJiq > m g also on-shell intermediate squark states will occur [see Fig. [24jb]. This 
case is treated in analogy to the previous example by regularizing the pole in S3. The 
integration over S3 is hidden in the angular integrations so that the technique must be 
described in some detail. To transform the angular integration to the integration over S3, 
we define a reference frame in which p 2 defines the z-axis: 



<i6>i sin 6»i / d9 2 -~ {A + Bcos9 1 + Csm9 1 cos0 2 ) 



(76) 
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The integrals over 9 2 are given by 
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for X = A + B cos 9\, Y = Csin9 1 , X < 0, X 2 > F 2 . Using this formalism, the 
total partonic cross-section can finally be rewritten as 
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The singularity in S3 for > can be regularized by introducing the non-zero squark 
width in the propagator. The separation of the lowest-order qq contribution and the 
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integration of the residual (principal-value) integral can be carried out in the same way 
as in the previous case. 

For identical particles in the final state, gluino or squark pairs, singularities can occur 
in both s 4 and S3. The existence of interference terms [see Fig. E5c], involving both types 



of singularities, demands some special care. Let us consider the gluino-pair production 
process as an example. In that case the interference terms can be treated by substituting 

1 V (-) + i7rS(s 4 ) (78) 



S4 — le \S4 

-> — 1 -=v( —± - ) - in6(a + b cos 9 1 ). (79) 

S3 + ze a + le + cos 61 \a + bcosUiJ 

The simultaneous singularities (s 4 — ie)~ 1 (s 3 + ie)" 1 do not require a subtraction proce- 
dure, since these configurations are kinematically suppressed. However, if nig > m g and 
s > (trig + nig) 2 in our example, the product of the two non-zero imaginary parts gives 
rise to a non-zero contribution to the cross-section. 



B.2 Hadronic Differential Cross- Sect ions 

For the hadroproduction of mixed squark-gluino final states, 

h^K x ) + h 2 (K 2 ) — q{ Pl ) + ~g(p 2 ) [+g(h)/q(k 3 )/q(k 3 )}, (80) 

we introduce the following hadronic variables: 

S=(K 1 + K 2 f (81) 
T x = (K 2 - p 2 ) 2 - m\ T g = (K 2 - p 2 ) 2 - m\ 



U x = {K x - p 2 ) 2 - m\ U g = {K x - p 2 f - m 2 . 



The double-differential hadronic cross-section is given by the integration of the parton 
cross-sections over the parton densities: 

(S,Tg,Ug,Q 2 ) = 



dTgdUg 



= E f ^ f dx 2Xl f^(x u Q 2 )x 2 f^(x 2 ,Q 2 ) tM^hlhl^l (82 ) 

i,j=g,q,q Jx l J x 2 M g dU g 

-T g - m? + m? _ -XiUg -m 2 ~+m 2 ~ 

Xl = S + U g %2 = Xl S + Tg ' 

Since the parton i has been assigned to the hadron hi, we can express the invariant parton 
variables in terms of the hadronic variables: 

S = X\X 2 S tg = X 2 Tg Ug = X\Ug. 
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The integration over x 2 can be transformed into the integration over S4, 

dxi j dx 2 = dxi / — — (83) 

JxZ Jx~ JO 1 lg 



x x (S + U g )+T g + ml-ml x 2 = - 



s 4 - X\U g -m 2 g + mj 
xiS + T g 



The integral over s 4 splits into contributions from soft and hard gluons. Both the soft and 
hard contributions contain terms that depend logarithmically on the cut-off parameter A. 
The soft terms of that type are mapped into the hard contributions, so that the scaling 
functions are independent of A and well defined in the limit A — > 0. 

The differential cross-section in transverse momentum and rapidity of the produced 
gluino can be obtained from the double-differential cross-section by the following trans- 
formation: 

dp t dy dT g dU g 

2 Tg Ug 2 tg Ug % ^ 1 I +g \ , s 

Pt = ~ m g = JL J^ L ~ m g y = 2 l0g [jj J ■ ( 85 ) 

[In general, transverse momentum and rapidity are defined for the particle carrying the 

momentum p 2 .] The transverse mass of the gluino is given by the quantity ^jp\ + m ~ . 

We define the rapidity spectrum by the sum of the contributions from positive and neg- 
ative valuesQ. For identical particles, the spectrum is defined for positive rapidity only, 
accounting in this way for the statistical factor 1/2 in those cases. Integrating over trans- 
verse momentum and rapidity, the total cross-section is reproduced: 

rpr ax (0) ry max (pt) d 2 a ry ma *(0) rpf ax (v) £ Q 

a(S, Q 2 ) = dp t dy —— = dy dp t -—- 

Jo Jo dptdy J Q J dp t dy , gg v 



Pt 



""•(v) = - - J {S + m 2 - m 2 ) 2 - 4m 2 ~S cosh 2 y 

2V o cosm/ v 



y max {p t ) = arccosh ' S + m ~9 



2JS(p 2 + ml 




In order to reproduce the total hadronic cross-section Eq. fl66|) , special care must be 
exercised in calculating the hadronic differential distributions for the resonance contribu- 
tions: 



13 Since squarks and antisquarks are not discriminated in the experimental analyses, we have not dis- 
tinguished them in the final states in the numerical discussion. For qq, qq, and gg final states we have 
therefore added the distributions with respect to both final-state particles. These y distributions are 
symmetric under y — > — y. Note that an extra factor 1/2 has to be inserted when using these (summed) 
distributions to obtain the total cross-section. 
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l/ s 4g singularity: For the subtraction it is convenient to write the double-differential 
hadronic cross-section in the form 



j2 rl rl 

dx! / dx 2 f^{x l ,Q 2 )f^{x 2 ,Q 2 



dTgdUg j T Jr/Xl 



J tg(s4) UZ 4 UL g 



= fdxj 1 dx 2 f^(x u Q 2 )f^(x 2 ,Q 2 ) r ds, 9 ^ 

Jt Jt/x 1 J0 S 4g 



J7) 



The integration boundaries are the ones defined in the previous subsection on the partonic 
processes (and r = (nig + nig) 2 / S) . The fact that the arguments of the 5-functions should 
have zeros inside the allowed integration region results in additional conditions on Xi, T g , 
and U g . As we have opted for a subtraction of on-shell intermediate states at the parton 
level, the subtracted double-differential hadronic distribution is given by 

''""^ fl dx 1 [ 1 dx 2 f*( Xl , Q 2 ) &(x 2 , Q 2 ) r ds, 9{S ' 9) - 9 S° 2 ] (88) 



(llgdUg J T Jt/x! JO b Ag^ m g L g 

after regularization by means of the non-zero gluino width. From this it is clear that upon 
integration over the hadronic variables T g and U g the subtraction for the total hadronic 
cross-section is reproduced. 

1/ss singularity: In this case we start off by applying Eq.([77|) and write 
d 2 a 



dTgdUg 



1 




dxi 


f 




Jt/xi 


r s t 






dS 3 g 


Jo 





39) 



> 4 ( s 3 9 ) f4(s4) J3i... / f \ / il 

xl <!,,„! dsj dtg a ° l \ + 5[Tg-^) 5[U g -^ 

ls^(s 3g ) Jt~{ Si ) ds 3g ds 4 dt g V x 2 J V xi 

Now the regularization by means of the non-zero squark width and subsequent subtraction 
of the resonance contribution can be performed in the usual way. 



Mixed l/s s and I/S4 singularities for identical final-state particles: After ex- 
changing the order of the integrations, the singularities in the propagators can easily be 
isolated. They can be treated as discussed in the subsection on partonic processes, giving 
rise to contributions from the two non-zero imaginary parts. 
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C Splitting Functions 



In this appendix we list all next-to-leading order SUSY-QCD splitting functions. The 
splitting functions involving massless partons [quarks and gluons] are known from the 
standard QCD evolution equations |35| , |36 , [25H . The splitting functions involving massive 



coloured SUSY particles are realized in final-state distributions at very high energies 



37| . We present all functions for the momentum fraction x in the restricted range < 



x < 1, excluding in this way the end-point singularities in P aa . The various end-point 
singularities can be derived from quark conservation, dx \P qq {x) + Pg L<1 (x) + Pq Rg (x)) = 

and Jq dx [P qq (x)+P qq (x)} = 0, and from momentum conservation, L dxxJ2 a Pab(x) = 0. 
[As usual, Cp = 4/3, iV = 3, and Tf = 1/2.] 



Quarks and gluons: 



Q 
9 

9 



q{+g) 
g(+q) 

q(+q) 

g(+g) 



Squarks and gluons: 

q -»• q(+g) 

q^ g (+<?) 
g -»• q (+q) 

Gluinos and gluons: 

g g(+g) 

9 -> g (+g) 
g^g (+g) 



P 9q (x) 
P q M 



Pgg( x ) 



c f 



1 + x 2 
1 — x 

l + (l-x) 



X 



T s [x 2 + (1 - x) 2 ] 
1 



2- N 



x(l — x) 



+ x(l-x)-2 



Pgq L ( x ) 
P9 l9 (x) 



P<i r q r ( x ) 



Pgq R ( x ) 



2C f 



x 



1 — x 
1 — X 

X 

P<i R g( x ) = T f x ( 1 - x ) 



2C r 



P rg {x) = N 
P gS (x) = N 



1 + x 2 
1 — X 

l + (l-x) 2 

X 



P 99 (x) = 2 • -N [x 2 + (1 - x) 2 ] 
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Squarks, gluinos and quarks: 



Q 
Q 

9 
9 
Q 
Q 



q(+~g) 
H+q) 

q(+q) 
q(+q L /q R ) 
H+~g) 
H+Ql/Qr) 



P M L 



Pq L q( X ) 



Pqq R ( x ) 
Pgq R ( x ) 



p~,A x ) 



Pq R q( x ) 



= C F 
= C F 

1 

= 2 TfX 
[l-x) 
1 



Cf x 



\c F {l-x) 



The splitting functions for the charge-conjugate transitions are identical. 



Using the notation of Section |3.2.3| , the end-point singularities in P aa are given by 
C F [2 log 5 + 1] 5(1 - x) for P qq and P m , and [2N\og5 + (3 /2] 5(1 - x) for P gg and P gg . 
Note that in Section |3.2.3| the end-point singularities for P gg and P qq have been derived 
for the case in which massive particles are decoupled. 
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